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Abstract 

A  coarse  measure  of  a  discrete  signal's  power 
distribution  vs.  frequency  is  its  "spectral  tilt", 
which  we  define  as  the  slope  of  least-squares 
linear  fit  to  the  log  power  spectrum 
(0  ^  CO  ^  71  rads/sec).  One  application  of  cal¬ 
culating  spectral  tilt  is  to  help  discruninate 
between  voiced  and  unvoiced/silent  tim<» 
segments  while  processing  sampled  speech 
signals.  When  linear  predictive  coding  ^PC) 
of  speech  is  used  to  model  its  short-time 
spectrum,  it  is  computationally  more  efficient 
to  calculate  spectral  tilt  directly  from  the  LPC 
model  parameters  instead  of  from  the  log 
power  spectrum.  In  this  paper  we  present 
methods  for  calculating  spectral  tilt  from 
cepstral  coefficients,  pole  values,  and  poly¬ 
nomial  coefficients.  These  methods  are  all 
based  on  calculating  spectral  tilt  as  a 
weighted  sum  of  cepstral  coefficients. 

Introduction 

Calculation  of  linear  spectral  tilt  has  been  of 
interest  to  scientists  in  various  disriplinftg  [i- 
2],  as  it  provides  a  coarse  measme  of  a  sig¬ 
nal's  power  distribution  vs.  frequency.  Given 
z  -domain  transfer  function  H{z)  we  define 
spectral  tilt  to  be  the  slope  of  line 

yK(£>)^m„ptGi-¥bopf  that  best  fits  the  log 
power  spectrum  lnj//(e^“)p  over  the  range 
0  ^  ^  71 ,  by  the  method  of  least-squares. 
That  is,  when  error  E  is  defined  as 

71  2 

E  =  \  {m©-f-h-ln|//(e-'‘“)|^  I  da, 

0  ' 
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E  is  minimum  when  b^b^pt  and  m^m^p,. 
In  the  continuous-frequency  case,  calculating 
spectral  tilt  using  this  formula  is  difficult  ex¬ 
cept  for  the  simplest  .  An  approxi¬ 

mation  to  the  spectral  tilt  is  usually  found 
from  M  samples  of  lnj//(e-'“)p,  taken  uni- 
fonnly  over  the  range  0  ^  ^  tt  .  Estimating 

spectral  tilt  in  this  manner  falls  within  the 
general  class  of  least-squares  polynomial  in¬ 
terpolation  [3],  and  estimation  error  is  re¬ 
duced  to  zero  as  M  -» <».  When  H{z)  is  in 
the  form  of  an  all-pole  model, 

G 

it  is  advantageous  to  find  from  the  N 
pole  values  directly,  without  first  having  to 
calculate  samples  of  hi|/f(e-'“)p.  This  not 
only  reduces  the  number  of  calculations  re¬ 
quired,  but  also  provides  an  exact  measure  of 
spectral  tilt  of  the  LPC  log  power  spectrum. 
The  rest  of  this  paper  deals  with  estimating 
i%pt  from  various  LPC  filter  parameters, 
such  as  the  pole  values  {Pi ,  P2  >  •  •  • .  Pv }  • 

Background 

Let  us  represent  the  spectral  tilt  of 
5(©)=ln|//(e-'“)p  using  functional  notation: 
m  =  Tilt {5 (©)} .  The  underlying  principle  to 
our  method  is  the  linearity  property  of  least- 
squares  spectral  tilt  calculation. 


That  is,  if  5(co)  =  + 0282  (co) ,  then 

=  Tilt{a,S,  ((o)+ ^2^2  (m)} 

=  aJ'ilt{S^  (co)}+  a^Tilt{S^  (co)}. 

Thus  when  5(a))  =  L5t(a)),  spectral  tilt  may 
be  found  as  the  sum  of  tilts  from  each  term 
5^((o): 

m  =  ][]r/7t{5^(Q))}  = 

k  k 

To  prove  that  the  linearity  property  holds  in 
this  case,  consider  the  following:  when  opti¬ 
mally  fitting  samples  of  S(co)  (0^a)^7t) 
with  a  first-order  polynomial  by  the  method 
of  least  squares,  b^pt  and  m^pt  may  be  found 
by  solving  the  following  system  of  normal 
equations  Ax  —  y: 
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Here  the  polynomial  being  defined  is 
—  fn^ptdi + bgpf ,  and  frequency  values 
Git  =  kK/{M-\)  (^'=o,l,■.•,M-l)  are  M 
uniformly-spaced  samples  of  co  over  the  in¬ 
terval  0  ^  CO  ^  71 .  Vector  y  may  be  ex¬ 
pressed  as  y  =  V^s,  where: 


'1  COq 

5(coi) 

1  a>i 

V  = 

1  0)2 

.1 

Since  matrix  multiplication  is  a  linear  opera¬ 
tion,  replacing  5  with  a  sum  of  terms  yields: 


N  n  /V 

=££>j. =£x,. 

/"I  {mI 


N 

I 
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jn^  +  m2  H - h  m^ 

Thus  we  have  proven  that  for  polynomial  fit¬ 
ting  of  sampled  points  by  the  method  of  least- 
squares,  the  sum  of  polynomial  coefficients 
calculated  to  optimally  fit  each  additive  com¬ 
ponent  of  s  equals  the  polynomial  coeffi¬ 
cients  calculated  to  optimtdly  fit  s .  The  true 
spectral  tilt  of  5(co)  is  approached  when  the 
number  of  sample  points  Af —>«>,  so  this 
property  is  also  valid  in  the  continuous  fre¬ 
quency  domain. 

Spectral  Tilt  from  Cepstral  Coefficients 
The  log  power  spectrum  5(a))  may  be  ex¬ 
pressed  as  a  summation  of  cosine  functions, 
weighted  by  the  cepstral  coefficients  [4]: 

5(a))  =  ln|//(e'“)f  =  21n|//(e^)| 

00  00 

=  cos(^'a))  =  T’5t(a)) 

i-l  IS1 

Thus  5(0)  is  the  sum  of  infinitely  many  co¬ 
sine  terms.  Using  the  linearity  property  dis¬ 
cussed  above,  we  may  express  spectral  tilt  in 
terms  of  the  cepstral  coefficients: 

00  eo 

C0S(^C0)} 

kml  i-l 

^  it-1,3.5.—  k 


Vector  X  may  be  found  after  a  single  matrix 
multiptication:  x=  A~V^s  =  Ds  (A~^  will 
always  exist  for  unique  sampled  values  of  co). 
Note  that  matrix  D  is  independent  of  5(co). 


Because  the  even-indexed  cepstral  coeffi¬ 
cients  are  the  weights  for  cosine  waves  hav¬ 
ing  an  integer  number  of  periods  in  the  range 
0  ^  CO  ^  7t  rads/sec,  there  is  zero  contribution 
to  spectral  tilt  from  these  terms. 


Spectral  Tilt  from  Linear  Predictor  Poles 
The  log  power  spectrum  of  an  N-th  order 
LPC  filler  may  also  be  expressed  in  terms  of 
its  N  pole  values  and  gain  constant  G : 


S{(o)  =  ]n\H{e^)f  =  In  Gn(l  - 

N 

=  lnG^-2£ln|l-p^e-^“| 


In  this  form  Sica)  is  the  sum  of  + 1  terms. 
The  spectral  tilt  of  the  term  InG^  equals 
zero,  leaving  the  following  expression  for 
spectral  tilt  of  S'(®)  in  terms  of  its  N  pole 
values: 

N 

m  =  r//r{S(©)}  *  -2£7’/7i{ln|l- 

We  were  not  able  to  solve  the  integral  and 
express  m  in  terms  of  the  pole  values 
{Pi,P2,--,Pn}  directly  from  this  equation. 
An  alternate  solution  method  exists,  how¬ 
ever:  to  first  calculate  cepstral  coefficients 
from  the  pole  values,  then  solve  for  spectral 
tilt  by  the  formula  derived  in  the  previous 
section.  The  cepstral  coefficients  corre¬ 
sponding  to  the  N-th  order  all-pole  log 
power  spectrum  may  be  found  as  follows  [5]: 


Combining  this  result  with  the  previous  ex¬ 
pression  for  spectral  tilt  in  terms  of  the  cep¬ 
stral  coefficients,  we  obtain  the  following  in¬ 
finite  summation: 


m  =  - 


-48  f 

^  *=1.3.5.- 


1  ^ 


n=l 


Fortunately,  due  to  the  p'^/k^  terms,  the 
solution  for  spectral  tilt  m  converges  rapidly 
with  increasing  k  (typically,  only  fifteen  or  so 
terms  yield  a  solution  to  full  TF.F.P.  double 


floating-point  precision).  Another  LPC  pole- 
based  method  is  to  pre-calculate  sjjectral  tilt 
for  various  combinations  of  pole  magnitude 
and  angle,  then  store  the  results  in  a  look-up 
table.  We  obtained  good  quantization  levels 
by  uniform  sampling  of  0  vs.  ln(l  -  r) ,  where 
p  =  re^ . 


Spectral  Tilt  from  Linear  Predictor 
Coefficients 

The  log  power  spectrum  of  an  A/-th  order 
LPC  filter  may  also  be  expressed  in  terms  of 
its  linear  predictor  coefficients  and  gain  con¬ 
stant  G : 


Sifo)  =  ln|/f(^-'®)|^  =  InG^  -In 


N 

1-  £at 
k=\ 


In  this  form  5(co)  may  not  be  expressed  as  a 
sum  of  terms,  where  each  term  depends  on  a 
coefficient  <Xj. .  The  approach  we  take  is  to 
take  advantage  of  a  recursive  solution  for  the 
cepstral  coefficients  fi*om  the  LPC  coeffi¬ 
cients  [6]: 

k-l 

r, 

^k-^  k-\ 

i  L  k>N. 

n=k-N 


The  cepstral  coefficients  are  then  used  to  cal¬ 
culate  spectral  tilt  from  the  formula  previ¬ 
ously  derived.  In  most  cases  this  method  will 
be  more  computationally  efficient  than  factor¬ 
ing  the  LPC  polynomial  to  obtain  the  pole 
values,  and  solving  for  spectral  tilt  using  the 
pole  equations. 

Experimental  Results 

We  have  applied  the  equations  derived  above 
to  measuring  the  spectral  tilt  of  sampled 
speech  data:  a  30  msec  harming  window  was 
stepped  across  the  speech  samples  with  0.1 
msec  steps,  the  windowed  data  zero-padded 
to  1024  points  prior  to  calculating  its  FFT. 


The  MATLAB  "polyfit"  routine  was  used  to 
estimate  the  spectral  tilt  of  513  samples  of 
ln|//(e^‘“)|^,  spaced  uniformly  between  0  and 
5  KHz.  A  180  ms  section  of  the  speech 
waveform  used  for  this  experiment  is  plotted 
as  Figure  1(a).  The  corresponding  spectral 
tilt,  calculated  directly  from  ln|i7(e-''“)p,  is 
plotted  as  the  solid  line  in  Figure  1(b).  (The 
slope  value  is  expressed  in  dB/kHz.)  This 
plot  shows  interesting  spectral  slope  features 
corresponding  to  silence,  onset  and  end  of 
voiced  speech,  sibilant  and  plosive  sounds. 


Figure  1(a).  A  180  ms  section  of  female  speech  se- 
lected  for  analysis. 

Figure  1(b).  Solid  line:  exact  spectral  tilt,  calculated 
from  the  log-magnitude  spectrum.  Dashed  line: 
spectral  tilt  estimated  from  the  poles  of  the  corre¬ 
sponding  UPC- 14  spectral  model. 


Next  we  calculated  LPC-14  coefficients  cor¬ 
responding  to  each  short-time  spectrum 
(autocorrelation  method),  factored  the  de¬ 
nominator  polynomial  to  obtain  pole  values 
for  each  frame,  and  applied  the  formula  pre¬ 
sented  above  that  calculated  LPC  spectral  tilt 
from  the  pole  values.  The  result  is  plotted  as 
the  dashed  line  in  Figure  1(b).  We  also  found 
that  lower  LPC  orders  provide  reasonable 
estimates  of  spectral  tilt,  as  shown  in  Fig.  2. 


Conclusion 

From  the  equations  presented  here  one  may 
calculate  spectral  tilt  directly  fi’om  various 
LPC  model  parameters,  avoiding  the  neces¬ 
sity  and  computational  expense  of  first  calcu¬ 
lating  the  log  power  spectrum.  We  rely  on 
the  linearity  property  of  least-squares  poly¬ 
nomial  data  fitting  to  simplify  our  expres¬ 
sions.  The  results  may  be  applied  to  effi¬ 
ciently  estimate  the  short-time  spectral  tilt  of 
speech  signals  to  help  discriminate  between 
various  speech  waveform  components. 


Figure  2.  Solid  line:  the  LPC-14  spectral  tilt  firom 

Figure  1(b).  Dashed  line:  LPC-2  spectral  tilt  for 

comparison. 
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